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ABSTRACT 

Recent radio observations of the small Galactic supernova remnant Gl •9+0-3 made at 
4.86 GHz with the VLA are presented, and compared with earlier observations at 1.49 GHz 
which have a comparable resolution (10 x 4 arcsec^). These show that the radio emission 
from this remnant has expanded significantly, by about 15 per cent over 23 years, with a cur- 
rent outer diameter of « 92 arcsec. This expansion confirms that Gl-9+0-3 is the youngest 
Galactic remnant yet identified, only about 150 years old at most. Recent, lower resolution, 
1.43-GHz observations are also discussed, and the integrated flux densities from these and the 
4.86-GHz observations are compared with earlier results. This shows that the integrated flux 
density of Gl •9+0-3 has been increasing recently. 
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1 INTRODUCTION 

The small Galactic 'shell' type supernova remnant (SNR) 
Gl-9+0-3 was identified by Green & Gull (1984), from Very Large 
Array (VLA) observations of a sample of Galactic radio sources. It 
was revealed as a small limb-brightened shell, brighter in the north, 
only just over 1 arcmin in diameter, and is the smallest catalogued 
Galactic SNR (Green 2004' ). Although its distance was not known 
at the time of its identification, it was clear that it must be a young 
SNR, due to its small physical size (e.g. even at distance as large 
as 17 kpc - i.e. twice the distance to the Galactic Centre ~ its phys- 
ical diameter would be ~ 6 pc, which is comparable to that of the 
remnant of Tycho's SN of AD1572). As such Gl-9+0-3 is one of 
few 'young but distant' SNRs, which are missing from current cat- 
alogues, due to selection effects (e.g. Green 2005), but perhaps also 
indicating a true deficit (The et al. 2006). Searches for these miss- 
ing young remnants have had limited success (see Green & Gull 
1984; Helfand et al. 1984; Green 1985, 1989; Sramek et al. 1992; 
Misanovic, Cram & Green 2002, see also Saikia et al. 2004). 

There are few published observations of GT9+0-3 in the lit- 
erature. VLA images at 4.9-GHz and 1.5-GHz (from observations 
made in 1985) are presented in Green & Gull (1984) and Green 
(2004) respectively. Gray (1994a) presents MOST observations at 
843 MHz, but these barely resolve the remnant. Nord et al. (2004) 
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' See also: http://www.mrao.cam.ac.uk/surveys/snrs/ for an up- 
dated version (from 2006 April). 



Table 1. VLA observations of Gl-9+0-3. 





1985 observations 


2008 observations 


date 


1985 April 16 


2008 Mai-ch 12 


VLA array 
frequency 


B 
1.49 GHz 


C 

4.86 GHz 


time on source 


25 min 


29 min 


primary calibrator 
assumed flux density 

secondary calibrator 
derived flux density 


3C286 

14.70 Jy 

B1829-106 

0.927 Jy 


3C286 

7.49 Jy 

J1751-251 

0.569 Jy 



and Yusef-Zadeh, Hewitt & Cotton (2004) present VLA observa- 
tions, at 330 MHz and 1.5 GHz respectively, but of limited qual- 
ity or resolution. Gl-9+0-3 has also been detected in X-rays, with 
ASCA as AX 11748.7-2709 (Sakano et al. 2002). 

Recently Chandra observations (Reynolds et al. 2008, here- 
after Paper I) have revealed that Gl-9+0-3 is one of the few 
shell SNRs that have line-free X-ray emission, dominated by syn- 
chrotron rather than thermal emission. The size of the X-ray emis- 
sion observed by Chandra was, unusually, significantly larger than 
the then published radio observations, which were interpreted in 
Paper I as requiring a large expansion rate - ~ 16 per cent over 23 
years - which implies that Gl -9+0-3 is a very young SNR. Here we 
present recent new radio observations made with the VLA, which in 
comparison with previous observations at a similar resolution, but 
a different frequency, confirm the large expansion rate, and hence 
youth of this remnant. We also discuss the temporal evolution of the 
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1,49 GHz observed 1985 April 16 



4,86 GHz observed 2008 March 12 




-27" 10' 00" - ^ 




Figure 1. VLA observations of Gl-9+0-3: (left) at 1.46 GHz from 1985 with contour levels at (-3, -2, -1, 1,2,3 . . . 10, 15,20. . . 95) x 0.30 mjy beam"', 
and (right) at 4.86 GHz from 2008 with contour levels at (—3,— 2,-1, 1,2,3... 10, 15,20... 95) x 0.17 mJy beam"'. The negative contours are dashed. Both 
images have a resolution of 10 x 4 arcsec^ at 5.5° W from N. The central crosses mark the centre position used for radial profiles (see text and Fig. 2). 



flux density of Gl-9+0-3, and conclude that it has been brightening 
over recent decades. 



2 OBSERVATIONS 

Gl-9+0-3 was recently observed with the VLA at 4.86 GHz in the 
C-array configuration. The resolution of these observations closely 
match those of VLA observations made in 1985 at 1.49 GHz in 
the B -array configuration (results from which are shown in Green 
2004). Details of these observations are given in Table 1. The new 
observations were processed using standard techniques in (classic) 
AIPS (see Bridle & Greisen 1994). Obviously corrupted data were 
flagged, the flux scale calibration was based on observations of 
3C286, and antenna-based amplitude and phase calibrations from 
observations of calibrator source(s) made nearby in time were ap- 
plied. The earlier 1.49-GHz observations were also reprocessed. 
Although the time spent observing Gl-9+0-3 was similar in each 
observation, the details of the Mv-coverage for the two observations 
are different, due to differences in LST coverage, missing antennas, 
and flagging of bad data. Consequently the synthesised beams are 
slightly different. For expansion studies comparison images were 
made at a matched resolution, see Fig. 1. These images show large 
expansion of the radio emission from Gl-9+0-3, which is discussed 
further in Section 3.1. The positional alignment of these images 
was checked from the apparent position of a compact source that is 
about 1.4 arcmin north of the centre of the remnant. The positions 
of this source were found to be in agreement to better than 0.4 arc- 
sec, which is small compared with the large changes seen between 
the 1985 and 2008 image of Gl-9+0-3. It should be noted, however, 
the earlier 1.49-GHz observations probably do not have sufficient 
short baselines to be able to fully image the largest scale emission 
from Gl-9+0-3. The VLA 'Observational Status Summary'^ lists 
the largest angular scale that 'can be imaged reasonably well' in 
a full synthesis as 2 arcmin for 1.49 GHz observations in the B- 
array, and notes that this value should be halved for single snap- 
shot observations. Thus, given the angular size of Gl-9+0-3 is over 
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1 arcmin, the largest scale structure may not be fully imaged in the 
1.49-GHz observations. For the 4.86-GHz C-array observations the 
largest angular scale well imaged with a full synthesis is 5 arcmin 
- since the C-array is not an exactly scaled version of B-array, but 
has better small baseline coverage - so even halving this for snap- 
shot observations, the full emission from Gl-9+0-3 should be well 
imaged. Indeed, from Fig.l, the 1.49-GHz shows extended negative 
artefacts, which are features of the limited small baseline coverage. 
The 4.86-GHz image is also better than the 1.49-GHz image, in 
terms of its sensitivity: the nominal noise on these images, from 
the r.m.s. away from bright emission, are about 0.038 and 0.012 
mJy beam"' at 1.49 and 4.86 GHz respectively. 

Given the problem with missing large-scale emission in the 
1.49 GHz observations, we also made some short 1.43-GHz ob- 
servations of Gl-9+0-3 in the same observation run as the 4.86- 
GHz observations, in order to investigate temporal variations in the 
flux density of Gl-9+0-3. These consisted of two short (Ri 4 min) 
scans of Gl-9+0-3, together with adjacent calibration observations 
of 3C286 to set the overall flux density scale (with an assumed 
flux density of 14.7 Jy), and of the nearby secondary calibrator 
J1751— 253 (for which a flux density of 1.18 Jy was derived). The 
results from these observations are discussed below in Section 3.2. 



3 DISCUSSION 

3.1 Expansion and Structure 

From Fig. 1 it is clear that there is a large expansion of the radio 
emission seen from Gl-9+0-3 between 1985 and 2008. To quan- 
tify this, radial profiles of the emission from each image were con- 
structed - see Fig. 2 - averaging over all angles, from a central po- 
sition which was adjusted to maximise the peak of the radial profile 
from the 2008 image. (This position is 17''48™45^4, -27° 10' 06", 
J2000.0, which close to the geometrical centre of the remnant.) A 
simple re-scaling of the 1985 radial profile indicates an average ex- 
pansion of 15 per cent between 1985 and 2008, i.e. 0.65 per cent 
yr" ' , is required to match the positions of the peaks of the radial 
profiles (Fig. 2). This expansion rate is biased towards that appro- 
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Chandra X-ray observed 2007 February/March 




Figure 2. Top: radial profiles of radio emission from Gl ■9+0-3, from 1985 
(dashed) and 2008 (solid). Bottom: scaled profile from 1985, expanded by 
12, 15 and 18 per cent (dotted, dashed, dotted), compared with the 2008 
profile. (The amplitudes of the profiles from 2008 have been multiplied by 
1.65.) 








Figure 3. The black contours, at (2.4,3.0,3.6... 7.8, 10.8, 13.8... 34.8) x 
10^** counts s^'beam"', are Chandra X-ray observations of Gl-9+0-3 (Pa- 
per I), with the same resolution and scale as the radio images shown in 
Fig. 1. (The background level away from the remnant is ps 1.8 x 10^** 
counts s^'beam^'.) The dotted blue contour is the lowest contour of the 
2008 4.86-GHz radio image shown in Fig. 1. 



priate to the brighter N and NE portions of the remnant, since the 
radial profiles are weighted by intensity. 

This result confirms that a significant expansion is needed 
to match the extent of the X-ray emission from the recent Chan- 
dra observations to older radio observations, as found in Paper I. 
Assuming free expansion, this implies an age of 150 years for 
Gl-9+0-3, which makes this the youngest known Galactic SNR 
(but also see Morris et al. 2006 and Renaud et al. 2006 for other 
possible very young SNRs which have been reported). As the rem- 
nant will have undergone some deceleration, this age is an upper 
limit. Very high X-ray absorption (Paper I) suggests that Gl-9+0-3 
is not much closer than the Galactic centre, while the very high 
expansion would indicate a very large velocity if it were much 
further away. Following Paper I, we adopt a nominal distance of 
8.5 kpc, which for an age of < 150 years and a outer radio diameter 
of 92 arcsec corresponds to a mean expansion speed of > 12,000 
km s^ . 

The interpretation of the radial profiles in simple terms of just 
a radial scaling is, however, complicated by a variety of effects. 
First, recall that the in'-coverage of these observations differ, par- 
ticularly - as noted above - that the 1.49-GHz observations may 
not well image the large scale emission from Gl-9+0-3. Second, 
the images are at different frequencies, so that any change in spec- 
tral index across the remnant would appear as a relative change in 
emission between the epochs (for example, a spectral index change 
of 0.05 would correspond to a relative intensity change of 6 per cent 
between 1.49 and 4.86 GHz). Third, there are real changes in 
the emission over the 23 years between the observations, so the 
changes seen are not just simple, self-similar expansion. In the case 
of Gl-9+0-3, we are not observing the expansion of distinct fea- 
tures, such as the radio knots in Gas A (see Anderson & Rudnick 
1995), or the well defined outer shock of Tycho's SNR (e.g. Tan & 
Gull 1985; Reynoso et al. 1997). From Fig. 1 there is evidence for 
changes in the shape and brightness of the radio emission between 
1985 and 2008. This is particularly noticeable in the NW, where 



there is a clear extension outside the brighter rim of emission from 
the remnant in the 2008 image (labelled 'NW in Fig. 1). There 
is some indication of extension in the 1985 image, but at a much 
lower level than in the 2008 image. The azimuth range of this outer 
extension corresponds closely to where the X-ray emission shows a 
bright 'ear' in this region (Paper I, see Fig. 3 above). In the E there 
is also an faint extension outside the main emission particularly in 
the 2008 image (labelled 'E'), but over a much smaller azimuth 
range than the outer extension in the NW. Also, the ridge line of 
the emission in the NE clearly differs between the epochs, being 
concave (w.r.t. to the centre of the remnant) in the 1985 image, but 
slightly convex in the 2008 image (labelled 'R'). If extrapolated, 
this apparently points towards the faint E extension, which corre- 
sponds to the northern edge of the X-ray emission outside the main 
rim in the E. Also, the peaks of the bright ridge of the radio emis- 
sion change between the two epochs (e.g. the peak in the NE in the 
1985 image splits into two peaks in the 2008 image). 

From Figs 1 and 3, there are striking differences between 
the overall structure of Gl-9+0-3 at radio and X-ray wavelengths, 
which may reflect differences in particle acceleration rates or effi- 
ciencies at different locations. Although the outline of the remnant 
is very nearly circular in the radio, the brightness of the radio emis- 
sion varies considerably in azimuth (with the bright NE rim be- 
ing about fifty times brighter than the weakest emission in the S). 
If a local density gradient were invoked to explain the large N-S 
difference in radio brightness, then it is not obvious why the rem- 
nant maintains such a circular outline. On the other hand, the X-ray 
emission of the remnant shows a bipolar structure, brightest in the 
E and W, which suggests the influence of a large scale magnetic 
field is important for the efficiency of accelerating the particles re- 
sponsible for the X-ray emission, but not for the those producing 
the radio emission. The brightest emission in the E and particularly 
in the NW corresponds to relatively faint radio emission. In the 
NW this is the outer extension seen in the 2008 radio image, which 
presumably corresponds to the region of the outer shock with the 
highest velocity. 
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Table 2. Radio flux densities for Gl-9+0-3. 



V 


S 


AS 


observation 


reference 


notes 


/MHz 


/Jy 


/Jy 


date(s) 






332 


2.84 


0.10 


1986-1989 


LaRosa et al. (2000) 




408 


1.18 


0.07 


1969-1971 


Clark & Crawford (1974) 


a 


843 


1.0 


0.05 


1985-1991 


Gray (1994a) 


b 


843 


0.986 


0.031 


1997-2007? 


Murphy et al. (2007) 




1400 


0.748 


0.038 


1993-1996 


Condon et al. (1998) 


c,d 


1425 


0.935 


0.047 


2008 


this work 


c 


2695 


0.440 


0.044 


1981-1984? 


Reich etal. (1984) 




4850 


0.236 


0.016 


1990 


Griffith et al. (1994) 


e 


4860 


0.437 


0.022 


2008 


this work 


c 


4875 


0.20 


0.05 


1974-1975 


Altenhoff etal. (1979) 


f 


5000 


0.20 


0.03 


<1975? 


CasweU et al. (1975) 


& 



Notes: 

(a) The listed flux density error is based on the quoted 5.4 per cent statistical 
flux scale (for 1 Jy sources), plus 2 per cent due to possible systematic 
declination variations, added in quadrature. 

(b) No flux density error is given, so a nominal value of 5 per cent is used; 
the observation dates are from Gray (1994b). 

(c) These flux densities were obtained from integrations of images, from 
within a polygon drawn around the remnant after removing a twisted plane 
fitted to pixels lying around the edge of the polygon (see Green 2007). The 
integrated flux density was found to vary by less than one per cent when 
a variety of polygons were used, but to be cautious we assume a statistical 
error of 5 per cent. 

(d)TheNVSS survey. Gl-9+0-3 is also listed in the NVSS source catalogue 
(Version 2.17): 58.8 x 33.3 arcsec^ in extent, with 0.7453 ±0.0248 Jy. 

(e) The flux density is from the Vizier version (see Ochsenbein, Bauer & 
Marcout 2000) of the PMN Tropical zone catalogue allowing for a general 
Gaussian width fit (see: Griffith & Wright 1993). The published flux density 
from a fixed widfli fit is 0.175 ± 0.014 Jy. 

(f) The flux densities of the sources in this survey are quoted to the nearest 
0.1 Jy, so a flux density error of 0.05 Jy is used. 

(g) The flux density error is from adding quoted errors of 10 per cent for 
uncorrected gain variations, 10 per cent for positional errors (which would 
lead to underestimates of the true value) and noise of 0.01 Jy, added in 
quadrature. 

3.2 Brightening 

Given that Gl-9+0-3 is such an apparently young SNR, then an ob- 
vious question is how is its radio flux density evolving with time? 
Table 2 hsts the available flux densities for G 1 -9+0-3 from the hter- 
ature, and from the new observations presented here. (There is also 
a flux density for Gl-9+0-3 at 365 MHz from the Texas survey, 
Douglas et al. 1996, but this has not been listed as it is not expected 
that it will represent the total flux of Gl-9+0-3, which is resolved 
by the instrument.) The available observations have been made with 
a variety of instruments and resolutions, and may not be on consis- 
tent flux density scales, making direct comparisons difficult. Nev- 
ertheless, there is evidence that the flux density of Gl-9+0-3 has 
been increasing over recent decades. 

(i) Gray (1994a) noted that his 843-MHz flux density was larger 
than expected from extrapolation between the earlier flux densities 
at 408 MHz and 5.0 GHz. Gray suggested that this discrepancy was 
due to the fact that the earlier 408-MHz and 5.0-GHz flux densities 
were derived assuming the source was unresolved by the ~ 3 ar- 
cmin beams of both sets of observations, and so the observed peak 
flux density would underestimate the integrated flux density. How- 
ever, even if the source were 1.2 arcmin in extent, as it was when 
identified by Green & Gull (1984), this effect would be less than 
about 10 per cent for an ~ 3 arcmin beam, which is not sufficient to 



explain this discrepancy. Also Gray noted that recent - i.e. presum- 
ably around 1994 - 4.85-GHz observations (from 'Haynes et al. in 
preparation', which we have not been able to identify in the litera- 
ture) give a flux density of 0.3 Jy for Gl-9+0-3, which is above the 
earlier 5.0-GHz flux density reported by Caswell, Haynes & Clark 
(1975). 

(ii) The 332-MHz flux density is considerably larger than the 
earlier value at 408 MHz (Clark & Crawford). (Note that LaRosa 
et al. quote a 1.4-GHz flux density from NVSS with an implau- 
sibly large error, i.e. 0.747 ± 0.250 Jy, and hence also give very 
large error in their derived 332 MHz to 1.4 GHz spectral index, i.e. 
0.93 ±0.23.) 

(iii) The integrated flux density of Gl-9+0-3 at 4.86 GHz from 
the observations described above is significantly larger than any of 
the values available in the literature at similar frequencies, even 
bearing in mind the large errors in some of the earlier flux densities 
at ~ 5 GHz, and the possibility of systematic differences in the flux 
density scales between the different observations. 

(iv) The 1.43-GHz flux density from our recent observations is 
about larger than the previously available value at 1.4 GHz from 
13 years earlier, by a factor of 1.25 ± 0.09, which implies a flux 
density increase of about ~ 2 per cent yr^ . 

Combining the integrated flux densities from the contemporaneous, 
new 1.45 and 4.86-GHz observations gives a spectral index, a, for 
the radio emission from Gl-9+0-3 (here defined in the sense that 
flux density scales with frequency as 5 o= v^") of 0.62. The error 
in this spectral index is ±0.06 using the nominal 5 per cent uncer- 
tainties in the individual flux densities, not including any possible 
systematic uncertainties in the relative flux density scales of the 
observations. However, is it difficult to reconcile this spectral index 
with all of the other observations in the literature. If Gl-9+0-3 is 
indeed brightening, then the spectral index between 332 MHz and 
1.4 GHz would be larger than the value of 0.93 obtained by LaRosa 
et al., as the 1.40-GHz was obtained after the 332-MHz data. Also, 
the two 843-MHz flux densities do not show obvious evidence for 
brightening. Excluding the 332-MHz flux density, it is possible to 
obtain a reasonable fit of a spectrum with a ~ 0.7 brightening at 
~ 2 per cent yr^' to the available flux densities for Gl-9+0-3. The 
model fits reported in Paper I assumed a flux density at 1 GHz of 
0.9 Jy. Increasing this to 1.2 Jy results in slightly steeper values of 
the radio-to-X-ray spectral index. However, the new values of both 
the spectral index and the rolloff frequency are still consistent with 
the earlier values at the 90 per cent confidence level. 

The fact that Gl-9+0-3 is brightening is not unexpected if it 
is indeed only a hundred years or so old. Radio observations of 
supernovae show brightening emission over time-scales of up to a 
year or so after the optical SN due to decreasing opacity. There- 
after these radio supernovae (RSN) show steady decline in emis- 
sion (see, for example, Weiler et al. 1986, 2002). Searches for radio 
emission from older RSN, up to about 80 years old (see Eck, Cowan 
& Branch 2002) mostly provide upper limits for the radio emis- 
sion, indicating that the luminosities of these older RSNe are less 
than that of the brightest known Galactic SNR. Thus, in order to 
evolve to match the observed properties of Galactic SNRs, it is ex- 
pected that radio emission from young SNRs will brighten, proba- 
bly on time-scales of about a century, when they have swept up suf- 
ficient interstellar medium (e.g. Gull 1973; Cowsik & Sarkar 1984). 
(SN1987A is brightening at radio wavelengths, as it interacts with 
ring of material very near the explosion; Manchester et al. 2002.) 
For a current angular size of 92 arcsec, and a 1-GHz flux density 
of 1.2 Jy, the surface brightness of Gl-9+0-3 is ^ 7.7 x 10"^" 
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W m^^Hz^^sr^'. This means that Gl-9+0-3 is brighter than the 
remnant of the SN of AD 1006, somewhat fainter than the remnants 
of Tycho's and Kepler's SNe of AD 1572 and 1604, but is consider- 
ably fainter than Cas A, which is otherwise the youngest (approxi- 
mately 300 hundred years old) and brightest known Galactic SNR 
(e.g. Green 2004). 

Predictions can be made for brightness increases based 
on simple assumptions about particle acceleration and possible 
magnetic-field amplification. For power-law behaviour of source 
radius with time, and of ejecta and ambient density with radius, 
power law time dependences can be found for the relativistic par- 
ticle energy density (assuming a constant efficiency, i.e. Ug o= pv,), 
magnetic field (assuming either flux freezing upstream or ampli- 
fication with a different efficiency), and therefore the synchrotron 
luminosity: 5v o= f'. Once substantial deceleration has begun (i.e. 
once a reverse shock has formed, within a few years of the explo- 
sion), almost all models decrease in luminosity with time. For con- 
stant ambient density the drop is slowest; for constant ambient mag- 
netic field as well, the luminosity can actually rise, roughly as t 
for the type la SN model of Chevalier (1982); all other cases drop 
with time. While these are simplistic estimates, they do demon- 
strate that a growing radio flux density requires special conditions. 

If Gl. 9-1-0.3 is now ~ 100 yr old, the observed flux increase 
by a factor of 1.25 over 13 years gives I = 1.6, far larger than the 
largest estimates above. While other explanations cannot be ruled 
out, the most natural explanation is that the efficiency with which 
shock energy goes into relativistic electrons and/or magnetic field 
must be increasing with time. 



4 CONCLUSIONS 

Comparison of VLA observations made in 1985 and 2008 confirm 
that Gl-9+0-3 has a large expansion rate of ~ 0.65 per cent yr^' . 
This means that this is a very young supernova remnant, at most 
only 150 years old, and hence the youngest known Galactic super- 
nova remnant. There is also evidence that the this remnant has been 
brightening over the last few decades at radio wavelengths, suggest- 
ing that the efficiency of particle acceleration and/or magnetic-field 
amplification has been increasing. Further high resolution, multi- 
epoch radio observations are required to study Gl-9+0-3's dynam- 
ics and structure in detail, and monitor the temporal evolution of its 
flux density. 
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